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ABSTRACT 

We present phase-referenced VLBI results on the radio continuum and the OH 18 cm megamaser 
emission from the Ultra-Luminous Infrared Galaxy, IRAS 17208-0014. The observations were carried 
out at 1599 MHz using the Very Long Baseline Array, the phased VLA, and the Green Bank Tele- 
scope. The highest resolution radio continuum results show several compact sources with brightness 
temperatures on the order of 10 6 K. These sources are more likely to be clustered supernova remnants 
and/or luminous radio supernovae. However, the agreement between the number of observed and 
expected compact sources above the 5 a level supports the possibility that each one of the compact 
sources could be dominated by a recently detonated luminous radio supernova. The continuum results 
suggest that there is no radio-loud AGN in the nuclear region of this galaxy. The OH 18 cm mega- 
maser emission in IRAS 17208-0014 is detected at various angular resolutions. It has an extent of 
170 x 110 pc, and is mostly localized in two regions separated by 61 pc. The structure and dynamics 
of the maser emission seem to be consistent with a clumpy, rotating, ring-like geometry, with the 
two dominant maser regions marking the tangential points of the proposed rotating-ring distribution. 
Assuming Keplerian motion for the rotating maser ring, the enclosed dynamical mass and the mass 
density within a radius of 30.5 pc, are about 3 x 10 7 (sin~ 2 i) M Q , and 281 (sin -2 i) Mq pc~ 3 , 
respectively. 

Subject headings: galaxies: individual (IRAS 17208-0014) — galaxies: starburst — radio continuum: 
galaxies — radio lines: galaxies — techniques: interferometric 



1. INTRODUCTION 

At bolometric luminosities greater than 10 11 L@, in- 
frared galaxies become the most numerous objects i n 
the local universe (z < 0.3) l|Sanders & Mirabel! ITflflfih . 
The trigger for the intense infrared emission appears to 
be the strong interaction or merger of molecular gas- 
rich spirals. Galaxies with the highest infrared lumi- 
nosities (£ir(8 — 1000 /im) > 10 12 L Q ), known as Ultra- 
Luminous Infrared Galaxies (ULIRGs) , appear to be ad- 
vanced merger systems powered by a nuclear starburst 
and/or an AGN. These may represent an i mportant stage 
in th e formation of quasi-stellar objects ijSanders et alJ 
1988). 

The extreme merging conditions in ULIRGs could, in 
some cases, induce nuclear OH 18 cm masing several 
orders of magnitude stronger than for typical Galactic 
masers, i.e., megamasers (10 6 times), or ev en gigamasers 
(10 9 ti mes) jBaanlll989t IBaan et al.llT992]) . Recent sur- 
veys bv lDarling fc Giovanellil l|2002) report an OH mega- 
maser (OHM) detection rate of up to 33% for ULIRGs 
at z > 0.1. 

Sensitive Very Long Baseline Interferometry (VLBI) 
observations provide the most direct means to study 
OHM and radio continuum emission in the innermost, 
dust-obscured, regions of merging systems. While radio 
continuum studies play a key role to directly image and 
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determine the nature of the nuclear power sources (AGN 
and/or starburst) in these galaxies, the OHM emission 
studies can enhance our understanding of the kinematics 
on parsec scales in their nuclear regions, and be used as 
probes of AGNs or intense nuclear starbursts. 

In the present paper, we report VLBI observations of 
the 18 cm radio continuum and OHM emission from the 
galaxy IRAS 17208-0014 at z = 0.04281. This object is 
one of the few strong OHM emitters at z < 0.1 that can 
be imaged at milliarcsecond (mas) resolution, and con- 
tribute to a better understanding of the nature, extent, 
and dynamics of extragalactic OHM emission. 

This galaxy, which is classified as a ULIRG, has an in- 
frared luminosi t y of L tr = 2.5 x 10 12 L Q , as defined in 
iGoldader et alJ |l995|) . We adopt a distance of 171 Mpc, 
assuming H Q = 75 km s -1 Mpc -1 . At this distance, 
1 mas corresponds to 0.76 pc. 

Optical images of IRAS 17 208-0014 at 6550 A show 
two tidal tails from a merger ijMelnick fc Mirabel Il99(l 
IMurphv et all 11996]) . Near-IR images show a very dis- 
turbed morphology and an extended, b ut single, nucleus, 
suggesting a very advanc ed merger ijZenner fc Lenzei] 
119931 IMurphv et all 119961) . Infrared observations sug- 
gest that this galaxy represents t he extreme of star burst- 
dominated sources of this type (iSoifer et alJl2000]) . and 
an observational proof that a collision of galaxies can 
lead to a mass distribut ion similar to elliptical galaxies 
ijZenner fc Lenzenlll993|) . 

The radio continuum emission of this galaxy at mil- 
liarcsecond (mas) res olution was firs t detec ted by the 
VLBI observations of iMomiian et all !|2( )()•"> ■ at 21 cm. 
These observations revealed that the radio emission is 
dominated by an extreme nuclear starburst, with no indi- 
cation of a compact high-brightness temperature source 
(i.e., a radio- loud AGN). 
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IRAS 17208- 0014 exhibits OHM activity at 1665, 1667, 
and 1720 MHz l|Martin et al]ll989fi . The strongest emis- 
sion is in the 1 667 MHz line, wi th a luminosity of 
L OH = 1054 L Q (jMartin et al.l|l9 89). Very Large Array 
(VLA) and single-dish Nangay observations showed the 
existence of two separate velocity structures with simi- 
lar w idths for both the 1667 and 1665 MHz OH main 
lines JMartin et al II 19891). However the 18 cm VLBI ob- 
servations of IDiamond et all l|1999|) did not detect the 
1665 MHz line, nor the continuum emission, and recov- 
ered only 20% of the 1667 MHz line. 

2. OBSERVATIONS AND DATA REDUCTION 

The present observations were carried out at 1599 MHz 
on 2003 May 25, using three NRAO 4 facilities: the 
Very Long Baseline Array (VLB A), the phased VLA, 
and the Green Bank Telescope. The bandwidth of the 
observations was 16 MHz, in each of the right- and left- 
hand circular polarizations, with the central frequency 
set appropriately for the 1667 MHz transition of the OH 
18 cm line at a heliocentric velocity cz = 12790 km s -1 
(jMartin et al.ll989tlMbmiian et alJ2003ft . The data were 
two-bit sampled, and correlated at the VLBA correlator 
in Socorro, New Mexico, with 512-point spectral resolu- 
tion per baseband channel, and 2 s correlator integration 
time. The total observing time was 10 hr. 

Data reduction and analysis were performed using 
the Astronomical Image Processing System (AIPS) 
and the Astronomical Information Processing System 
(AIPS++). 

Along with the target source IRAS 17208-0014, the 
compact source J1730+0024 was observed as a phase ref- 
erence. A cycle consisting of 180 s on the target source 
and 60 s on the reference source was adapted. The source 
J1743-0350 was used for amplitude and bandpass cali- 
bration. 

After applying a priori flagging and manually excis- 
ing integrations affected by interference, we performed 
amplitude calibration using the measurements of the an- 
tenna gain and system temperature of each station, and 
bandpass calibration. The flux density calibration of the 
data set is believed to b e accurate at the level of 5% 
HWrohel TJlvestadll2005T). 

A continuum data set was generated by averaging the 
OH-emission-free spectral channels. In the continuum 
data set, the phase reference source J1730+0024, was 
self-calibrated and imaged in an iterative cycle. The 
self-calibration solutions for J 1730+0024 were applied to 
both the continuum and line data of the target source. 

The continuum and line data sets of IRAS 17208-0014 
were then deconvolved and imaged at various spatial 
resolutions, with a grid weighting intermediate between 
pure natural and pure uniform (ROBUST = in AIPS 
task IMAGR). The continuum images were obtained us- 
ing the multi-resolution CLEAN algorithm in IMAGR, 
which has been shown to be superior compared to the 
standard CLEAN or Maximum Entropy Method in imag- 
ing a combinatio n of very weak exten ded and compact 
emission regions l|Momiian eTaT1l27)0! . 

3. RESULTS 

4 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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Fig. 1. — Continuum image of the central region in IRAS 17208- 
0014 at 1599 MHz. The restoring beam size is 17.6 X 6.1 mas at 
position angle —4°. The peak flux density is 217 £tjy beam -1 , and 
the contour levels are at —3, 3, 4, 5 and 6 times the rms noise level, 
which is 31 /^Jy beam -1 . The gray-scale range is indicated by the 
wedge at the right side of the image. The reference position (0, 
0) is a(J2000) = 17 h 23 m 21?9648, <5(J2000) = -00°17W819. For 
comparison purposes, this origi n has been chosen to be the same 
as for the continuum figures of lMomiian et aT H2003T) . 



Fig. 1 is the continuum image of the central region 
in IRAS 17208-0014 at 1599 MHz and 17.6 x 6.1 mas 
(13.4 x 4.6 pc, P.A. =—4°) resolution, and represents 
th e bright emissio n regio n to the west of center in Fig. 2 
of lMomiian et alJ l)2003|) . This region is populated by 
10 compact sources with peak flux densities > ha 
(155 ,uJy beam -1 ), but < 7a (217 /^Jy beam -1 ). Most of 
these sources were reproduced with pure natural weight- 
ing (ROBUST=5), but not with pure uniform weighting 
(ROBUST=-5). This is largely due to the much higher 
noise in the uniformly weighted image, which is a factor 
of 6 worse than natural weighting, and almost a factor of 
4 worse than the intermediate weighting (ROBUST=0). 
Table 1 lists the Gaussian fitting parameters of the com- 
pact continuum sources seen in Fig. 1 derived using the 
AIPS task JMFIT. The positions (col. [2]) are relative 
to the first source in the list. Column (3) lists the sur- 
face brightnesses of these sources, and column (4) their 
integrated flux densities. Column (5) gives the decon- 
volved sizes of the Gaussian components at FWHM, and 
column (6) gives the position angles. For several of these 
components, the nominal deconvolution sizes are listed 
as given by JMFIT. However, for some components, no 
nominal sizes were obtainable for the minor axes of the 
fitted Gaussians, hence, the maximum sizes of the mi- 
nor axes are reported as upper limits. Moreover, because 
source number 10 is unresolved, we are only able to quote 
an upper limit for its overall size. The corresponding 
brightness temperatures of these compact sources are on 
the order of 10 6 K. 

Fig. 2 is a total intensity (moment 0) image of the OH 
18 cm main lines (both 1665 and 1667 MHz) at 25 mas 
resolution. The image was obtained by blanking areas 
where the flux density of the OH spectral- line image cube 
is below 4 er, or 3.6 mJy beam -1 . Fig. 2 also shows the 
Hanning-smoothed spectra of the OH 18 cm main lines 
averaged over various regions of the OHM emission. 

OH emission from IRAS 17208-0014 is spread over a 
broad range of velocities. The strong, narrow feature 
seen in all profiles of Fig. 2 (other than for Region 4) al- 
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most certainly arises from the 1667 MHz transition. This 
conclusion seems very likely for two reasons. Firstly, this 
relatively narrow feature lies close to the systemic veloc- 
ity of the galaxy if the feature arises from the 1667 MHz 
transition (see Section 2). Secondly, this narrow feature 
is always accompanied by a weaker emission feature off- 
set by about 370 km s -1 , the velocity equivalent to the 
difference in frequency between the 1665 and 1667 MHz 
lines at the redshift of IRAS 17208-0014. 

The overall extent of the megamser emission is 170 x 
110 pc on the plane of the sky, and seems to be domi- 
nated by two regions separated by 80 mas (61 pc) at po- 
sition angle 83°. We designate these as the eastern and 
western components (Regions 1 and 5 in Fig. 2). These 
two regions seem to be connected by weaker maser emis- 
sion regions, interrupted in the middle by what seems 
to be a region devoid of maser emission. In addition to 
the two main maser regions, a weak extension toward 
the south, and an even weaker one toward the north are 
detected. The measured velocity difference between the 
eastern and the western regions is 137 ± 3 km s -1 . 

Table 2 summarizes the physical characteristics of the 
dominant OH peaks seen in the spectra of Fig. 2 for 
both the 1665 and the 1667 MHz transitions. The line 
properties were obtained by fitting Gaussians to the non- 
Hanning-smoothed OH spectra obtained in the regions 
specified in Fig. 2 using the tool IMAGEPROFILEFIT- 
TER in AIPS++. Column (1) lists the region number 
as designated in Fig. 2. The velocities (col. [2]) refer 
to the center velocities of the 1667 MHz line peaks in 
each region. Column (3) lists their full widths at half 
maximum, and column (4) is their peak flux density val- 
ues. Columns (5)-(7), are similar to columns (2)-(4), 
but for the 1665 MHz line peaks. Column (8) is the peak 
flux density ratio of the 1667 MHz line to the 1665 MHz 
line for each region. For the spectrum of region Re- 
gion 6, which shows two distinct emission peaks for the 
1667 MHz line, two Gaussian components are listed. No 
Gaussian fitting parameters are reported for the spec- 
trum seen in Region 4, because of its low signal-to-noise 
ratio and the apparent blending of the two mainlines. 

Fig. 3 is a Hanning-smoothed spectrum of the OH 
18 cm main lines integrated over the whole OH emis- 
sion region seen in Fig. 2. The velocity scale is refer- 
enced to the 1667 MHz line. The weak emission be- 
tween 130 00 and 13400 km s" 1 is from the 1665 MHz 
main line ijMartin et al.lll989fl . The arrow in this fig- 
ure points to the expected velocity of the 1665 MHz line 
that corresponds to the main peak of the 1667 MHz line. 
While the shape of the spectral profile seen in Fig. 3 
is similar to that ob tained at Nancay and with the VLA 
ijMartin et al.ll989|) . the VLBI array recovers only ~80% 
of the single-dish and VLA flux densities. 

Fig. 4 displays gray-scale total intensity (moment 0) 
images of the 1667 MHz main line at three different an- 
gular resolutions with continuum contours superposed. 
For comparison p urposes, the same orig in is used for the 
coordinates as bv lMomiian et all (|2003f) for their contin- 
uum images. Fig. 4a is at 50 mas resolution, Fig. 46 at 
25 mas resolution, and Fig. 4c at the full resolution of 
the VLBI array, which is 17.6 x 6.1 mas resolution (P. 
A.= —4°). These images were obtained by restricting 
the spectral channel range to that of the 1667 MHz main 
line, and blanking the OH image cubes at the 4 a level, 




Fig. 2. — Gray-scale and contour total intensity (moment 0) 
image of the OH 18 cm main lines, and Hanning-smoothed OH 
spectra obtained at various locations of the OHM emission. The 
restoring beam size is 25 X 25 mas. The contour levels are at 0.3, 
0.6, 0.9, . . . , 1.8 Jy beam -1 km s — 1 , and the gray scale range is 
indicated by the step wedge at the right-hand side of the image 
in units of Jy beam -1 km s — 1 . The reference position (0, 0) is 
a(J2000) = 17 h 23 m 21?9533, <5(J2000) = -00° 17' 00!' 983. The rms 
noise level of the total intensity is 0.13 Jy beam -1 km s — 1 . The 
velocity scale of the spectra is referenced to the 1667 MHz line, 
and the spectral featur es between 13000 a nd 13400 km s — 1 are 
from the 1665 MHz line IMartin et alJ19891) . The effective velocity 
resolution of the spectra is 12.3 km s -1 
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Fig. 3. — Hanning-smoothed OH 18 cm spectrum integrated 
over the complete OH emission region as seen in our VLBI obser- 
vations. The velocity scale is referenced to the 1667 MHz line. The 
effective velocity resolution is 12.3 km s — 1 . The arrow points to 
the expected velocity of the 1665 MHz line that corresponds to the 
main peak of the 1667 MHz line. At the redshift of IRAS 17208- 
0014, the offset between these two lines is 368 km s — . 

where a is 1.8, 0.9, and 0.5 mJy beam -1 , for the 50 mas, 
25 mas, and the full resolution images, respectively. In 
all these images, the OH emission from the eastern re- 
gion appears to be stronger than that from the west. 
The lower resolution images show the eastern region to 
be associated with stronger continuum emission than the 
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Fig. 4. — Gray-scale total intensity (moment 0) images of the 
1667 MHz OH 18 cm main line, with 1599 MHz continuum contours 
superposed, at various angular resolutions, (a) At 50 mas resolu- 
tion: Contour levels are at 4, 8, 12, 16, 20 X la (89 fijy beam -1 ). 
(6) At 25 mas resolution: Contour levels are at 4, 6, 8, 10 X 
1<t (57 £ijy beam -1 ). (c) At the full resolution of the array 
(17.6 X 6.1 mas, P.A.=-4°): Contour levels are at 4, 5, 6 X la 
(31 /ijy beam -1 ). The reference position (0, 0) is the same as in 
Fig. 1. The gray scale range in all these images is indicated by 
the step wedge at the right-hand side of each image in units of 
Jy beam -1 km s — 1 . 

western. Moreover, the continuum emission seems to 
have a larger extent than that of the OHM, and its peak is 
located north-east of the megamaser regions. At the full 
resolution, we find no correlation between the brightest 
continuum and maser emission spots. The derived lower 
limit on the maser amplification factor is ~ 130. 

Fig. 5 is the velocity field (moment 1) of the 1667 MHz 
main line at 50 mas resolution. A general velocity gra- 
dient is seen in the east-west direction, however, the ve- 
locity contours seem to be distorted toward the western 
region. Moreover, the eastern region itself shows a north- 
south velocity gradient. 

4. DISCUSSION 
4.1. The Radio Continuum 
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Fig. 5.— The OH 1667 MHz velocity field (moment 1) at 50 mas 
resolution. Contours are from 12700 to 12850 km s -1 in steps 
of 25 km s -1 . The gray-scale range is between 12670 and 12910 
km s — 1 . The reference position (0, 0) is the same as in Fig. 2. 

The radio continuum emission from the ULIRG 
IRAS 17208-0014 at 1362 MH z has been discussed in de- 
tail by IMomiian et al.l l|2003f) for various angular resolu- 
tions. Our observations at 1599 MHz have imaged the ra- 
dio continuum emission from IRAS 17208-0014 at a con- 
siderably higher angular resolution tha n previously re- 
portcd, which was 36 x 33 mas, P.A. 25° (|Momiian et alJ 
2003). 

The continuum image at the full resolution of our VLBI 
array (i.e., 17.6 x 6.1 mas; Fig. 1) shows several compact 
sources with flux densities greater than 5 times the rms 
noise level of 155 fiJy beam -1 . We find no dominant 
compact radio source with a very high brightness tem- 
perature. The brightness temperature of these sources 
are on the order of 10 6 K, indicating that the emission is 
non-thermal. 

I ■ the discussion in Momjian. et al. (2003) and 

references therein, these compact sources are likely to 
be clustered radio supernovae (RSNe) and/or super- 
nova remnants (SNRs), with the likelihood that each of 
these sources could be mainly powered by a recently- 
detonated, individual, RSN. Using the po st-maximum 
light curve of luminous RSNe derived by iSmith et alJ 
(1998), we find that individual RSNe will lie above our 
5 a detection level for 3.6 yr. Gi ven the supernova rat e 
of 4±1 yr- 1 in IRAS 17208-0014 ijMomiian et al.ll200^) . 
we expect to see about 14 ± 4 individual RSNe. This is 
in good agreement with the 10 sources revealed in our 
observations. 

We defer further discussion of the continuum emission 
to a subsequent paper, in which we hope to report lower 
angular resolution results at 1599 MHz with matching- 
resolution lower-frequency VLBI continuum observations 
at 330 and 610 MHz, and higher frequency MERLIN 
continuum observations at 5 GHz. 

4.2. The OHM Emission 

The spectral-line observations reported here reveal the 
details of the OHM emission from IRAS 17208-0014. 
The maser emission is complex, and seems to be localized 
in two distinct areas, the eastern and western compo- 
nents, separated by 61 pc (Fig. 2). Weaker maser emit- 
ting regions can be seen toward the north and the south. 

Our VLBI observations recover ~ 80% of the single 
dish flux density of the O HM emission in IRAS 17208- 
0014 ijMartin et al.lll989|) . This suggests that extended 
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OH emission is resolved by even the shortest baseline in 
our VLBI array, namely the phased- VLA to Pie Town 
baseline (52 km), and indicates that such emission re- 
gions have angular extents of at least 0"37. This is well 
within the ~ 5" extent of t he galaxy as seen in multi- 
wavelength IR observations l|Soifer et al.ll2006T) . 

The overall nature, geometrical distribution, and kine- 
matics of the megamaser emission in IRAS 17208-0014 
show striking simila rities to those of another megamser 
galax y, III Zw 35 ijPihlstrom et alJ 120011: IParra et all 
2005). These authors explain the OH emission in 
III Zw 35 by invoking the presence of an inclined clumpy 
ring structure oriented in a north-south direction. A sim- 
ilar ring model, though with an east-west orientation, 
could describe the OH emission in IRAS 17208-0014. In 
such a model, the two strongest emission peaks, the east- 
ern and western components (Fig. 2), would mark direc- 
tions tangential to the ring, and the apparent empty re- 
gion in the middle would coincide with the center of the 
ring. Assuming a circular ring structure, and excluding 
the southern and northern regions seen in Fig. 2, the ap- 
parent extend of the OHM emission on the plane of the 
sky suggests an inclination angle of ~ 55°. Numerical 
simulations and modeling, however, are needed to derive 
a more reliable estimate for the inclination angle. 

In addition to the eastern and western components, 
the OHM emission in IRAS 17208-0014 shows a promi- 
nent extension to the south at various angular resolu- 
tions. This structure does not seem to be part of the pro- 
posed ring structure, but lies perpendicular to it. Con- 
sidering its position, width, and velocity (Fig. 2, Region 
7 and Table 2), this feature could be associated with 
an inflow or an outflow. Molecular outflows and in- 
flows have been previously reported in several megamaser 
galaxies, sometimes with veloci t ies up to 800 km s" 1 
ijBaan. Haschick. fc Henkelll98flt IPihlstrom et al.l l2005V 
Simulations carried out by iMihos fc Hernauisn (1996) 
show that when gas-rich galaxies merge, strong gravi- 
tational torques drive the gas to the center of the rem- 
nant as a strong inflow and fuel an intense but very 
short-lived starburst. Considering that IRAS 17208- 
0014 is an advan ced merger system powered by an ex- 
treme starburst ( |Zenngr fc Lenzenlll993t iMurphv et al.1 
H996tlMomiian et all2003l) . the southern extension could 
be the result of an inflow that is fueling the star bust. Al- 
ternatively, the southern extension could be associated 
with an outflow. Extreme merging conditions, and con- 
sequently powerful nuclear activity (AGN and/or star- 
burst), could induce superwinds in the nuclear region, 
and cause material to be blown away. However, our 
data cannot uniquely establish that the southern exten- 
sion is an inflow or an outflow. Considering that an- 
other, weaker, extension is seen to the north (Region 2 
of Fig. 2), it might be the case that both the north- 
ern and southern extensions are part of another rotating 
structure, or they may represent different maser spots in 
a turbulent velocity field. 

The larger extent of the continuum emission relative to 
that of the OHM emission in IRAS 17208-0014 (Fig. 4a- 
b) indicates that the star formation in the nuclear re- 
gion of this galaxy is not confined to just around or 
within the proposed ring structure, suggesting that the 
OH population inversion occurs only over a small re- 
gion. This is likely because the inversion of the OH 



main lines depends critically on the local IR spectrum 
(See Parra et al. 2005, and references therein). Moreover, 
these two images show that the OHM and peak contin- 
uum intensities do not coincide. This can be explained 
by free-free absorption , as was suggested for III Zw 35 
IPihls trom et al] 120011) . If the optical depth is signifi- 
cant at 1.6 GHz (r > 1), then the emission measure, 
EM = nil > 6 x 10 6 pc cm" 6 . Assuming I ~ 100 mas 
(76 pc), then the density is n e > 281 cm -3 , which is 
consis tent with that expected f or a dense starburst nu- 
cleus (Sanders & Mirabel 1996). Our subsequent paper 
on multi-frequency VLBA and MERLIN continuum ob- 
servations will address the likelihood of free-free absorp- 
tion in IRAS 17208-0014. Another possibility is that 
this system is an incomplete merger, where the nuclei of 
the progenitor galaxies have not yet completely merged. 
The OH ring traces one of the original nuclei, and the 
star formation in the center is asymmetrically distributed 
around this nucleus. 

At the full resolution of the array, the lack of posi- 
tional coincidence between the continuum and megamser 
spots in IRAS 17208-0014 (Fig. 4c) is also found in 
other OH M galaxies, e.g . Arp 220, III Zw 35 and 
Mrk 273 iLonsdale et all [19981 IPihlstrom et~aTI 120011 
iKlockner fc Baanll2004|) ." The derived lower limits on the 
maser amplification fa c tors i n IR AS 17208-0014 are up 
to 130. lLonsdale et al.l l| 19981) and lDiamond et all l|1999h 
estimate that the amplification factors of the megamser 
emission in Arp 220 and III Zw 35 are up to 800 and 
500, respectively. Furthermore, these authors explain 
the lack of positional correlation between the continuum 
and the OHM spots by satur ated compact masing clouds. 
However. IParra et al.l |2005). and based on their detailed 
modeling of the OHM emission in III Zw 35, propose in- 
stead that the masers could be unsaturated. The bright- 
est maser emission could occur in regions where the con- 
tinuum emission is either very weak, or below the current 
detection threshold, because of the exponential effect of 
path le ngth on unsaturat ed masers occurring in clumpy 
media (|Parra et alJl2005[) . 

The measured 1667 to 1665 MHz flux density ratio 
at various regions of the OHM emission varies between 
3.6 ± 0.4 and 8.4 ± 0.5 (Table 2). Similar, high 1667 
to 1665 MHz ratios are reported in other OHM galax- 
ies, such as Mrk 273 and III Zw 35 ijYates et al J 1200(1 
Pihlstrom et alJl2001[) . These high ratios clearly indi- 
cate that the OH molecular levels are not in local ther- 
modynamic equilib rium (LTE), wh ich would yield a ra- 
tio of 1.8 l|Henkel fc WilsorJll990)h Moreover, the line 
ratios of the eastern and western OHM regions (Re- 
gions 1 and 5; Fig. 2 and Table 2) in IRAS 17208-0014, 
which are the tangents of the proposed ring structure, are 
very similar to those o f the tangent points in III Zw 35 
IjPihlstrom et alJl200H) . Based on the ve ry high 1667 to 
1665 M Hz line ratios observed in Arp 220. lLonsdale et al] 
( 1998) have suggested the existence of two maser emis- 
sion phases. One is compact, characterized by saturated, 
high gain, masers, and the other diffuse, characterized by 
unsatu rated, low gain, masers. However. lPihlstrom et alJ 
l)2001[) have proposed instead a mechanism for OHM 
emission based on a single phase of OH masing small 
clouds (~1 pc), and IParra et"aTI |2005) have argued 
that the observed differences in the line ratios between 
compact and diffuse phases could be a natural prop- 
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erty of one phase characterized by clumpy unsaturated 
m asers, in agre ement with the classical megamaser model 
of iBaanl (|j989). Furthermore, our results suggest that 
higher line ratios seem to be correlated with slightly 
narrower 1667 MHz line widths. Regions with line ra- 
tios greater than 6 show velocity widths between 73 and 
85 km s _1 , while regions with lower line ratios display 
velocity widths between 95 and 128 km s _1 . Therefore, 
higher 1667 to 1665 MHz line ratios can be explained 
by the exponential amplification of the background radi- 
ation by overlappin g unsaturated maser clouds in both 
space and velocity IjParra et al.N2005t ILockett fc Elitzurl 

The exponential amplification that causes stronger 
1667 MHz lines, and consequently higher 1667 to 1665 
MHz line ratios, could also explain the narrower mea- 
sured widths of the 1667 MHz lines compared to those 
at 1665 MHz (Table 2). While the 1667 MHz transition 
exhibits a range of strong and weak features in a given 
region (Figure 2), the 1665 MHz counterparts will have 
a smaller range of amplitudes, and appear as a wider, 
blended, spectral feature. 

The similarities between the OHM emission in 
IRAS 17208-0014 and III Zw 35, and consequently 
the proposed ring model, are further supported by the 
observed velocity field in IRAS 17208-0014 (Fig. 5). 
This field shows a strong gradient in the eastern region 
in two perpendicular directions, east- west, and north- 
south, and a distorted field in the western region. A 
similar velocity field h as been observed in III Zw 35 
l)Pihlstr5m et al.ll2001|) . where the northern masing re- 
gion has a two directional velocity gradient, and the 
southern region has a distorted velocity field. 

To explain the observed velocity fiel d by the ring mode l 
for the OHM emission in III Zw 35. iParra et~aLl (|2005f) 
have proposed that masing clouds, in addition to their 
rotational velocity, have a comparable outflow velocity 
parallel to the ring axis and directed away from the ring 
midplanc. In the terminology of these authors, a cloud 
near a tangent region that lies above the ring midplanc 
will have a projected outflow velocity that is blueshifted, 
while a cloud below the ring midplane will have a red- 
shifted projected outflow velocity. Such a mechanism 
could produce an apparent velocity difference of nearly 
twice the projected outflow velocity, and would be per- 
pendicular to the rotational velocity of the ring. These 
authors further note that at the tangent point where ma- 
terial is moving away from the observer, i.e., the north- 
ern tangent point in III Zw 35, or the eastern tangent 
point in IRAS 17208-0014, the cloud outflow mechanism 
gives a velocity gradient in the same direction as that 
caused by the rotation mechanism, thus reinforcing the 
gradient. In contrast, the tangent point where material 
is moving toward the observer, i.e., the southern tan- 
gent point in III Zw 35, or the western tangent point in 
IRAS 17208-0014, the two mechanisms have opposite di- 
rections and will partly cancel. This could result in the 
distorted velocity field seen in the southern component of 
III Zw 35, or in the western compo nent of IRAS 17208- 
0014. However, IParra et"aTI l)2005f) have suggested that 
alternatively the difference may mainly be due to the sta- 
tistical nature of the clumpy maser model used in their 
Monte-Carlo simulations. 

While various regions of the OHM emission show one 



dominant spectral peak in the 1667 MHz line, Region 6 
(Fig. 2 and Table 2) displays a double spectral feature, 
suggesting the existence of a velocity substructure in this 
region. Considering that bright OHM emission spots are 
the result of spatial and velocit y alignment of a num- 
ber of unsaturated maser clouds IjParra et al.ll2005lL the 
observed double-peak spectrum indicates that there are 
two distinct aligned groups of masing clouds within this 
region. 

Assuming a rotational velocity of ~ 69 (sin i) km s _1 , 
the velocity gradient between the eastern and west- 
ern masing peaks in IRAS 17208-0014, i.e., the 
tangents of the proposed ring distribution, is ~ 
2.3 (sin i) km s _1 pc -1 . For a radius of 30.5 pc, 
we derive an enclosed dynamical mass of about 
3 x 10 7 (sin~ 2 i) Mq. A greater mass value is estimated 
on la rger scales from mol ecular (CO) and atomic (H i ) 
gas (pownes fc Solomon! Il998t IMomiian etaH I2003|) . 
Within a radius of 30.5 pc, the measured mass density 
from the OH observations is 281 (sin -2 i) M Q pc~ 3 . 

The large scale velocity field attributed to the rotat- 
ing ring structure is comparable or smaller than the 
1667 MHz line widths obtained at various regions of the 
ring (Table 2). Very similar rotational velocity and ve- 
locity widths are repor ted in III Zw 35 ijPihlstrom et alJ 
l200lHParra et al.l2~0 05). These velocity widths, although 
very wide compared to galactic OH masers, seem to be 
typical for OHM sourc es. In their VLBI study of Arp 220, 
ILonsdale et alJ ljl998j) suggest that much of the broad 
megamaser line wi dth is intrin si c to t he compact maser 
spots. Moreover, IParra et al.l 1)2005(1 explain that the 
observed line widths cannot arise from a population of 
narrow-width clouds, but from the spectral blending of 
masing clouds with internal velocity widths 3> 1 km s _1 . 
In their detailed modeling of the OHM emission from 
III Zw 35, these authors report that in order to fit the 
observed velocity width of the compact maser spots, their 
model required maser clouds with internal velocity dis- 
persion of 20 km s _1 . The high velocity dispersion of the 
OHM clouds indicates a highly turbulent medium, and 
suggests that the c louds are not gravitationally confined 
IjParra et al.ll2005^ . 

5. CONCLUSIONS 

We have presented the results of sensitive, phase- 
referenced VLBI observations of the 18 cm continuum 
and OH megamaser emission in the central 0.4 kpc of 
the ULIRG IRAS 17208-0014. 

The continuum emission from this galaxy was imaged 
at a higher resolution than has been previously reported. 
Several compact sources with flux densities greater than 
5cr were detected. These sources have brightness tem- 
peratures on the order of 10 6 K, and are more likely to 
be clustered RSNe and SNRs. However, we cannot rule 
out the possibility that each of the compact sources is 
mainly powered by an individual bright RSN nested in 
a region that contains faded SNRs. This is supported 
by the agreement in the number of the observed and 
expected continuum sources above the 5 a level. The 
continuum results do not reveal a radio-loud AGN in the 
nuclear region of this galaxy. 

The OHM results show the detection of both the 1665 
and 1667 MHz mainlines in IRAS 17208-0014. The OHM 
emission in this ULIRG has an extent of 170x 110 pc, and 
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TABLE 1 

COMPACT SOURCES IN IRAS 17208-0014 



Gaussian Component Parameters 



Source 


Relative Position 3 


Peak b 


Total 


Deconvolved Sizc c 


P.A. 




(mas) 


(/ijy beam -1 ) 


(pJy) 


(pc) 


(°) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1 


0, 


177 ±31 


320 ± 81 


9.3 x 5.0 


5 


2 


13W, 32S 


196 ± 30 


604 ± 119 


21.3 x 5.4 


12 


3 


24W, US 


164 ± 32 


241 ± 71 


9.8 x < 7.1 




4 


57W, 60S 


168 ± 31 


299 ± 81 


7.6 x 4.4 


38 


5 


65W, 31S 


185 ± 31 


295 ± 75 


10.7 x 3.0 


10 


6 


81W, 36S 


167 ±31 


307 ± 83 


10.1 x < 8.9 




7 


131W, 29S 


166 ± 32 


243 ± 71 


12.3 x < 2.3 




8 


146W, 55S 


164 ± 32 


210 ± 66 


6.8 x 2.4 


7 


9 


185W, 18S 


184 ± 31 


319 ± 79 


7.1 x 5.4 


177 


10 


200W, 61S 


167 ± 33 


119 ±45 


< 4.9 





a The reference position (0,0) is a(J2000) = 17 h 23 m 21?96158, <5(J2000) = -00°17'00;'8904. 
b Higher than bo = 155 /ijy beam -1 . 
c At half maximum. 



is mostly localized in two regions separated by 80 mas 
(61 pc). The overall structure and dynamics of the maser 
emission seem to be consistent with a clumpy, rotating, 
ring-like geometry, with the two dominant maser regions 
marking the tangential points of the proposed ring struc- 
ture. However, our OH results reveal another component 
extended to the south. Considering the ongoing intense 
starburst in this advanced merger system, the southern 
extension could to be associated with a gaseous inflow or 
outflow. Assuming Keplerian motion for the observed 
OH ring structure, the enclosed dynamical mass and 
the mass density within a radius of 30.5 pc, are about 
3 x 10 7 (sin~ 2 z) M Q , and 281 (sin -2 i) M Q pc" 3 , re- 
spectively. 

Our results show that the starburst activity, as seen 
in the radio continuum images, has a larger extent than 
the OHM emission. Moreover, the low-resolution images 
show that at larger scales, the OHM emission peaks do 
not coincide with the peak of the continuum emission. 
This could be the results of free- free absorption and/or 



possibly an indication that the merger is incomplete. At 
the full angular resolution, the brightest continuum and 
OHM spots are not correlated. The derived lower limits 
on the maser amplification are up to 130. 

The observed line-to-continuum ratios, 1667 to 
1665 MHz line ratios, and the 1667 MHz line widths sug- 
gest that the strong OHM emission in various regions can 
be explained by unsaturated masing clouds overlapping 
in both space and velocity. 
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TABLE 2 

GAUSSIAN PARAMETERS OF THE OH MEGAMASER LINES IN THE REGIONS OF FIGURE 2 



1667 MHz 1665 MHz 



Region Velocity AVfwhm ^1667 Velocity AVfwhm i*i665 ^1667/^1665 

(km s _1 ) (km s -1 ) (mJy beam -1 ) (km s _1 ) (km s _1 ) (mjy beam -1 ) 

(1) (2) (3) (4) (5) (6) (7) (8) 



1 12864 ±2 73 ±4 12.3 ± 0.5 13193 ± 18 234 ± 40 2.0 ± 0.2 6.2 ± 0.5 

2 12861 ±8 97 ±13 4.8 ± 0.4 13217 ± 17 261 ± 44 1.3 ±0.2 3.7 ± 0.5 

3 12727 ± 2 85 ±4 11.3 ±0.3 13095 ± 13 270 ± 34 1.9 ±0.2 6.0 ± 0.4 

5 12727 ±2 77 ±4 13.4 ± 0.5 13103 ± 15 263 ± 39 1.6 ±0.2 8.4 ± 0.5 

6 12768 ± 2 95 ±5 8.0 ± 0.3 13140 ± 8 198 ± 18 2.2 ± 0.2 3.6 ±0.4 

12875 ± 2 56 ±4 6.1 ±0.4 

7 12724 ±2 128 ±5 6.5 ± 0.2 13114 ± 11 139 ± 27 1.2 ±0.2 5.4 ± 0.3 

8 12715 ± 3 82 ±5 9.9 ± 0.4 13083 ± 10 95 ± 24 1.2 ±0.3 8.3 ± 0.5 



